Datasets:

• AUG: Augmentation Data
• MAT: Matched Data
• AUG*: cleaned up AUG
• MAT*: Filtered MAT
• AUG* + MAT*

[a,b] = findgroups(osrm_data_sp.arac_id);

color_1 = rand(length(b),3);

dir_chosen = 1;

h = figure();

ax1 = subplot(2,1,1);

scatter(ax1, osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3,osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).speed ,5,color_1(a(osrm_data_sp.dir == dir_chosen)),'filled');

hold on

scatter(ax1, aug_data(aug_data.dir == dir_chosen,:).distance_from_start * 1e-3,aug_data(aug_data.dir == dir_chosen,:).speed ,5,"red",'filled','MarkerFaceAlpha',0.3);

title(ax1,['MAT + AUG --> dir = ', num2str(dir_chosen), ' speed color: MAT(arac-id) red:AUG'])

xlabel(ax1,'st km from kızılay')

ylabel(ax1,'speed km/hr')

yLimits = get(ax1,'YLim');

line(ax1,1e-3*[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

% text(ax1, 1e-3*markers.distance_from_start(markers.dir == dir_chosen), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

min_x = min(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

max_x = max(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

xlim(ax1,[min_x , max_x]);

ax1 = subplot(2,1,2);

scatter(ax1,osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3,osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).speed ,5,color_1(a(osrm_data_sp.dir == dir_chosen)),'filled');

hold on

scatter(ax1, aug_data(aug_data.dir == dir_chosen,:).distance_from_start * 1e-3,aug_data(aug_data.dir == dir_chosen,:).speed ,5,"red",'filled','MarkerFaceAlpha',0.3);

ylim(ax1,[-20 400])

yLimits = get(ax1,'YLim');

line(ax1,1e-3*[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3*markers.distance_from_start(markers.dir == dir_chosen), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

title(ax1,'y lim [-20 400]')

xlabel(ax1,'st km from kızılay')

ylabel(ax1,'speed km/hr')

min_x = min(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

max_x = max(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

xlim(ax1,[min_x , max_x]);

set(h,'Units','normalized','Position',[0 0 1 0.5]);

dir_chosen = 2;

h = figure();

ax1 = subplot(2,1,1);

scatter(ax1, osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3,osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).speed ,5,color_1(a(osrm_data_sp.dir == dir_chosen)),'filled');

hold on

scatter(ax1, aug_data(aug_data.dir == dir_chosen,:).distance_from_start * 1e-3,aug_data(aug_data.dir == dir_chosen,:).speed ,5,"red",'filled','MarkerFaceAlpha',0.3);

ax1.XDir = "reverse";

yLimits = get(ax1,'YLim');

line(ax1,1e-3*[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

% text(ax1, 1e-3*markers.distance_from_start(markers.dir == dir_chosen), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

title(ax1,['MAT + AUG <-- dir = ', num2str(dir_chosen), ' speed color:MAT(arac-id) red: AUG']);

xlabel(ax1,'st km from kızılay')

ylabel(ax1,'speed km/hr')

min_x = min(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

max_x = max(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

ax1.XTick = int64(linspace(min_x,max_x,9));

ax1.XTickLabel = num2cell(int64(linspace(min_x,max_x,9)));

xlim(ax1,[min_x , max_x]);

set(ax1, 'XTickLabel', flip(ax1.XTickLabel));

ax1 = subplot(2,1,2);

scatter(ax1,osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3,osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).speed ,5,color_1(a(osrm_data_sp.dir == dir_chosen)),'filled');

hold on

scatter(ax1, aug_data(aug_data.dir == dir_chosen,:).distance_from_start * 1e-3,aug_data(aug_data.dir == dir_chosen,:).speed ,5,"red",'filled','MarkerFaceAlpha',0.3);

ylim(ax1,[-20 400])

ax1.XDir = "reverse";

yLimits = get(ax1,'YLim');

line(ax1,1e-3*[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3*markers.distance_from_start(markers.dir == dir_chosen), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title(ax1,'y lim [-20 400]')

xlabel(ax1,'st km from kızılay')

ylabel(ax1,'speed km/hr')

min_x = min(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

max_x = max(osrm_data_sp(osrm_data_sp.dir == dir_chosen,:).distance_from_start * 1e-3);

xlim(ax1,[min_x , max_x]);

ax1.XTick = int64(linspace(min_x,max_x,9));

ax1.XTickLabel = num2cell(int64(linspace(min_x,max_x,9)));

set(ax1, 'XTickLabel', flip(ax1.XTickLabel));

% set(ax1, 'XTickLabel', flip(ax1.XTickLabel));

[a,b] = findgroups(osrm_data_sp.arac_id);

color_1 = rand(length(b),3);

dir_chosen = 1;

h = figure();

ax1 = subplot(2,1,1);

index = osrm_data_sp.filter & (osrm_data_sp.dir == dir_chosen);

index2 = aug_data.filter & (aug_data.dir == dir_chosen);

scatter(ax1, osrm_data_sp(index,:).distance_from_start * 1e-3,osrm_data_sp(index,:).speed ,5,color_1(a(index)),'filled');

hold on

scatter(ax1, aug_data(index2,:).distance_from_start * 1e-3,aug_data(index2,:).speed ,5,"red",'filled','MarkerFaceAlpha',0.3);

title([' MAT* + AUG* --> dir = ', num2str(dir_chosen), ' speed color:MAT(arac-id) red: AUG']);

xlabel('st km from kızılay')

ylabel('speed km/hr')

yLimits = get(ax1,'YLim');

line(ax1,1e-3*[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3*markers.distance_from_start(markers.dir == dir_chosen), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

min_x = min(osrm_data_sp(index,:).distance_from_start * 1e-3);

max_x = max(osrm_data_sp(index,:).distance_from_start * 1e-3);

xlim(ax1,[min_x , max_x]);

dir_chosen = 2;

ax1 = subplot(2,1,2);

index = osrm_data_sp.filter & (osrm_data_sp.dir == dir_chosen);

index2 = aug_data.filter & (aug_data.dir == dir_chosen);

scatter(ax1, osrm_data_sp(index,:).distance_from_start * 1e-3,osrm_data_sp(index,:).speed ,10,color_1(a(index)),'filled');

hold on

scatter(ax1, aug_data(index2,:).distance_from_start * 1e-3,aug_data(index2,:).speed ,5,"red",'filled','MarkerFaceAlpha',0.3);

hold off

title([' MAT* + AUG* <-- dir = ', num2str(dir_chosen), ' speed color:MAT(arac-id) red: AUG']);

xlabel(ax1, 'st km from kızılay')

ylabel(ax1,'speed km/hr')

min_x = min(osrm_data_sp(index,:).distance_from_start * 1e-3);

max_x = max(osrm_data_sp(index,:).distance_from_start * 1e-3);

xlim(ax1,[min_x , max_x]);

ax1.XDir = 'reverse';

yLimits = get(ax1,'YLim');

line(ax1,1e-3*[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3*markers.distance_from_start(markers.dir == dir_chosen), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

xlim(ax1,[min_x , max_x]);

ax1.XTick = int64(linspace(min_x,max_x,9));

ax1.XTickLabel = num2cell(int64(linspace(min_x,max_x,9)));

set(ax1, 'XTickLabel', flip(ax1.XTickLabel));

set(h,'Units','normalized','Position',[0 0 1 0.5]);

h = figure();

scatter(osrm_data_sp.space_distance,osrm_data_sp.speed,10,'filled')

xlabel('space difference b/w points (meter)')

ylabel('corresponding speed km/hr')

grid on

title(' MAT delta-x vs speed for both dir-s')

scatter(osrm_data_sp.time_distance * 60 * 60,osrm_data_sp.speed,10,'filled')

xlabel('time difference b/w points - seconds-')

ylabel('corresponding speed km/hr')

grid on

title(' MAT delta-x vs speed for both dir-s')

scatter(osrm_data_sp(osrm_data_sp.filter,:).space_distance,osrm_data_sp(osrm_data_sp.filter,:).speed,10,'filled')

xlabel('space difference b/w points (meter)')

ylabel('corresponding speed km/hr')

grid on

title(' MAT* delta-x vs speed for both dir-s')

scatter(osrm_data_sp(osrm_data_sp.filter,:).time_distance * 60 * 60,osrm_data_sp(osrm_data_sp.filter,:).speed,10,'filled')

xlabel(' time difference b/w points - seconds-')

ylabel('corresponding speed km/hr')

grid on

title(' MAT* delta-x vs speed for both dir-s')

scatter3(osrm_data_sp(osrm_data_sp.filter,:).time_distance * 60 * 60,osrm_data_sp(osrm_data_sp.filter,:).space_distance,osrm_data_sp(osrm_data_sp.filter,:).speed,10,'filled')

xlabel('time diff')

ylabel('space diff')

zlabel('speed')

grid on

title('MAT* dx-dt-v')

scatter(osrm_data_sp.space_distance(osrm_data_sp.speed > 200),osrm_data_sp.time_distance(osrm_data_sp.speed > 200) * 60 * 60,10,'filled')

h = figure();

dir_chosen = 1;

index = (osrm_data_sp.dir == dir_chosen) & osrm_data_sp.filter;

index2 = (aug_data.dir == dir_chosen) & aug_data.filter;

ax1 = subplot(2,1,1);

boxchart(ax1,double(osrm_data_sp.segment_id(index)) ,osrm_data_sp.speed(index),"MarkerStyle",'.',"WhiskerLineColor","blue",'LineWidth',1.5);

hold on

boxchart(ax1,double(aug_data.segment_id(index2)) ,aug_data.speed(index2),"MarkerStyle",'.',"WhiskerLineColor","red",'LineWidth',0.5);

legend(ax1,'AutoUpdate','off')

legend(ax1,'MAT*',"AUG*")

title(ax1,'MAT* + AUG* --> direction = 1 8:00-10:00 segment id vs speed boxplot')

xlabel(ax1,'segment id (kızılay - > polatlı)')

xlim(ax1,[min(double(osrm_data_sp.segment_id(index))) - 1, max(double(osrm_data_sp.segment_id(index))) + 1])

ylabel(ax1,'speed km/h')

% set(ax1,'YLim',[0 200])

% axis tight

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0.8*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

dir_chosen = 2;

index = osrm_data_sp.dir == dir_chosen & osrm_data_sp.filter;

index2 = (aug_data.dir == dir_chosen) & aug_data.filter;

ax2 = subplot(2,1,2);

boxchart(ax2,double(osrm_data_sp.segment_id(index)) ,osrm_data_sp.speed(index),"MarkerStyle",'.',"WhiskerLineColor","blue",'LineWidth',1.5);

hold on

boxchart(ax2,double(aug_data.segment_id(index2)) ,aug_data.speed(index2),"MarkerStyle",'.',"WhiskerLineColor","red",'LineWidth',0.5);

legend(ax1,'AutoUpdate','off')

legend(ax1,'MAT*',"AUG*")

min_x = min(double(osrm_data_sp.segment_id(index)));

max_x = max(double(osrm_data_sp.segment_id(index)));

ax2.XLim = [ min_x - 1 , max_x + 1];

% ax2.XTick = int64(linspace(min_x,max_x,10));

ax2.XDir = "reverse";

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title(ax2,'MAT* +AUG* <-- direction = 2 8:00-10:00 segment id vs speed boxplot')

xlabel(ax2,'segment id (kızılay -> polatlı)')

ylabel(ax2,'speed km/h')

% ylim(ax2,'auto')

% axis tight

yLimits = get(ax2,'YLim');

line(ax2,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax2, double(markers.segment_id(markers.dir == dir_chosen)), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

h=figure();

ax1 = subplot(2,1,1);

dir_chosen = 1;

index = aug_data.dir == dir_chosen;

index2 = osrm_data_sp.dir == dir_chosen;

[GC,GR] = groupcounts(aug_data.segment_id(index));

[GC2,GR2] = groupcounts(osrm_data_sp.segment_id(index2));

bar(ax1,GR2,GC2,'FaceAlpha',0.3,'EdgeAlpha',0.7);

hold on

bar(ax1,GR,GC,'FaceAlpha',0.3,'EdgeAlpha',0.7);

title(ax1,'AUG +MAT --> dir = 1 speed data count vs segment id')

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

xlabel('segment id')

legend(ax1,"MAT","AUG")

ax1 = subplot(2,1,2);

% [~,b2]=ismember(segments.segment_id,GR2);

% % segments.segment_id(nonzeros(b2));

% % segments.calc_length(nonzeros(b2));

% A = [GC2 , segments.calc_length(nonzeros(b2)) * 1e-3];

% A1L = [0; cumsum(A(:,2))];

% for k1 = 1:size(A,1)

% patch(ax1,[0 1 1 0]*A(k1,2)+A1L(k1),[0 0 1 1]*A(k1,1), "blue", 'LineWidth',0.01,'FaceAlpha',0.3,'EdgeAlpha',0.7);

% end

% hold on

[a,~]=ismember(segments.segment_id,GR);

A = [GC , segments.calc_length(a) * 1e-3];

% A1L = [0; cumsum(A(:,2))];

A1L = segments.distance_from_start(a) * 1e-3;

for k1 = 1:size(A,1)

patch(ax1,[0 1 1 0]*A(k1,2)+A1L(k1),[0 0 1 1]*A(k1,1), "red", 'LineWidth',0.01,'FaceAlpha',0.3,'EdgeAlpha',0.3);

end

[a2,~]=ismember(segments.segment_id,GR2);

A2 = [GC2 , segments.calc_length(a2) * 1e-3];

A1L2 = segments.distance_from_start(a2) * 1e-3;

for k1 = 1:size(A2,1)

patch(ax1,[0 1 1 0]*A2(k1,2)+A1L2(k1),[0 0 1 1]*A2(k1,1), "blue", 'LineWidth',0.01,'FaceAlpha',0.3,'EdgeAlpha',0.3);

end

axis(ax1,[min(A1L) max(A1L) 0 max(ylim)])

yLimits = get(ax1,'YLim');

line(ax1,1e-3 * [markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3 * double(markers.distance_from_start(markers.dir == dir_chosen)), 0.75*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title(ax1,'AUG --> dir = 1 speed data counts vs segments - bar width = segment length(km)')

xlabel(ax1,'km from kızılay')

ylabel(ax1,'aug data count')

h=figure();

ax1 = subplot(2,1,1);

dir_chosen = 2;

index = aug_data.dir == dir_chosen;

index2 = osrm_data_sp.dir == dir_chosen;

[GC,GR] = groupcounts(aug_data.segment_id(index));

[GC2,GR2] = groupcounts(osrm_data_sp.segment_id(index2));

bar(ax1,GR2,GC2,'FaceAlpha',0.3,'EdgeAlpha',0.7);

hold on

bar(ax1,GR,GC,'FaceAlpha',0.3,'EdgeAlpha',0.7);

title(ax1,'AUG +MAT --> dir = 2 speed data count vs segment id')

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

xlabel('segment id')

legend(ax1,"MAT","AUG")

ax1.XDir = "reverse";

ax1 = subplot(2,1,2);

[a,~]=ismember(segments.segment_id,GR);

A = [GC , segments.calc_length(a) * 1e-3];

% A1L = [0; cumsum(A(:,2))];

A1L = segments.distance_from_start(a) * 1e-3;

for k1 = 1:size(A,1)

patch(ax1,[0 1 1 0]*A(k1,2)+A1L(k1),[0 0 1 1]*A(k1,1), "red", 'LineWidth',0.01,'FaceAlpha',0.3,'EdgeAlpha',0.3);

end

[a2,~]=ismember(segments.segment_id,GR2);

A2 = [GC2 , segments.calc_length(a2) * 1e-3];

A1L2 = segments.distance_from_start(a2) * 1e-3;

for k1 = 1:size(A2,1)

patch(ax1,[0 1 1 0]*A2(k1,2)+A1L2(k1),[0 0 1 1]*A2(k1,1), "blue", 'LineWidth',0.01,'FaceAlpha',0.3,'EdgeAlpha',0.3);

end

min_x = min(A1L);

max_x = max(A1L);

axis([min_x max_x 0 max(ylim)])

set(h,'Units','normalized','Position',[0 0 1 0.5]);

ax1.XDir = "reverse";

yLimits = get(ax1,'YLim');

line(ax1,1e-3 * [markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3 * double(markers.distance_from_start(markers.dir == dir_chosen)), 0.75*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

xlim(ax1,[min_x , max_x]);

ax1.XTick = int64(linspace(min_x,max_x,9));

ax1.XTickLabel = num2cell(int64(linspace(min_x,max_x,9)));

set(ax1, 'XTickLabel', flip(ax1.XTickLabel));

title('AUG <-- dir = 2 speed data counts vs segments - bar width = segment length(km)')

xlabel('st km from kızılay')

ylabel('aug data count')

dir_chosen = 1;

index = (aug_data.dir == dir_chosen) & (~isnan(aug_data.speed)) & aug_data.filter;

[G,b] = findgroups(aug_data.segment_id(index));

h=figure();

ax1 = subplot(2,1,1);

plot(ax1, b, splitapply(@mean,aug_data.speed(index),G),'-+')

get_ave = @(x) hours(sum(x)) / length(x);

[~,y] = ismember(b,segments.segment_id);

hold on

plot(ax1, b, 1e-3 * segments.calc_length(y) ./ splitapply(get_ave, aug_data.time_diff(index),G),'-*')

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title({'AUG* direction = 1 8:00-10:00 segment id vs AVERAGE speed',...

"kantara bak"})

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

ax1.XLim = [ min(b) - 1 , max(b) + 1];

xlabel('segment id')

ylabel('speed km/h')

legend('mean v over segment','v from travel time mean')

dir_chosen = 2;

index = (aug_data.dir == dir_chosen) & (~isnan(aug_data.speed)) & aug_data.filter;

[G,b] = findgroups(aug_data.segment_id(index));

ax1 = subplot(2,1,2);

plot(ax1, b, splitapply(@mean,aug_data.speed(index),G),'-+')

get_ave = @(x) hours(sum(x)) / length(x);

[~,y] = ismember(b,segments.segment_id);

hold on

plot(ax1, b, 1e-3 * segments.calc_length(y) ./ splitapply(get_ave, aug_data.time_diff(index),G),'-*')

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title(ax1,' AUG* <-- direction = 2 8:00-10:00 segment id vs AVERAGE speed')

ax1.XDir = "reverse";

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

xlabel(ax1,'segment id')

ylabel(ax1,'speed km/h')

ax1.XLim = [ min(b) - 1 , max(b) + 1];

legend('mean v over segment','v from travel time mean')

dir_chosen = 1;

index = (aug_data.dir == dir_chosen) & (~isnan(aug_data.speed)) & aug_data.filter;

[G,b] = findgroups(aug_data.segment_id(index));

b2 = ismember(segments.segment_id,b);

h=figure();

ax1 = subplot(2,1,1);

plot(ax1, 1e-3*segments.distance_from_start(b2), splitapply(@mean,aug_data.speed(index),G),'-+')

get_ave = @(x) hours(sum(x)) / length(x);

[~,y] = ismember(b,segments.segment_id);

hold on

plot(ax1, 1e-3*segments.distance_from_start(b2), 1e-3 * segments.calc_length(y) ./ splitapply(get_ave, aug_data.time_diff(index),G),'-*')

title({'AUG* direction = 1 8:00-10:00 AVERAGE speed vs st km',...

"kantara bak"})

yLimits = get(ax1,'YLim');

line(ax1,[1e-3*markers.distance_from_start(markers.dir == dir_chosen),1e-3*markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3*double(markers.distance_from_start(markers.dir == dir_chosen)), 0*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

min_x = min(1e-3*segments.distance_from_start(b2));

max_x = max(1e-3*segments.distance_from_start(b2));

xlim(ax1,[min_x , max_x]);

% xlabel('segment id')

ylabel('speed km/h')

xlabel('st km from kızılay')

legend('mean v over segment','v from travel time mean')

%

ax1 = subplot(2,1,2);

dir_chosen = 2;

index = (aug_data.dir == dir_chosen) & (~isnan(aug_data.speed)) & aug_data.filter;

[G,b] = findgroups(aug_data.segment_id(index));

b2 = ismember(segments.segment_id,b);

plot(ax1, 1e-3*segments.distance_from_start(b2), splitapply(@mean,aug_data.speed(index),G),'-+')

get_ave = @(x) hours(sum(x)) / length(x);

[~,y] = ismember(b,segments.segment_id);

hold on

plot(ax1, 1e-3*segments.distance_from_start(b2), 1e-3 * segments.calc_length(y) ./ splitapply(get_ave, aug_data.time_diff(index),G),'-*')

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title({'AUG* direction = 2 8:00-10:00 AVERAGE speed vs st km'})

yLimits = get(ax1,'YLim');

line(ax1,[1e-3*markers.distance_from_start(markers.dir == dir_chosen),1e-3*markers.distance_from_start(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, 1e-3*double(markers.distance_from_start(markers.dir == dir_chosen)), 0*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

% ax1.XLim = [ 1e-3*min(segments.distance_from_start(b2)) - 1 , 1e-3*max(segments.distance_from_start(b2)) + 1];

ax1.XDir = "reverse";

min_x = min(1e-3*segments.distance_from_start(b2));

max_x = max(1e-3*segments.distance_from_start(b2));

xlim(ax1,[min_x , max_x]);

% ax1.XTick = int64(linspace(min_x,max_x,10));

% ax1.XTickLabel = num2cell(ax1.XTick);

set(ax1, 'XTickLabel', "");

% xlabel('segment id')

ylabel('speed km/h')

xlabel('st km from kızılay')

legend('mean v over segment','v from travel time mean')

dir_chosen = 1;

index = aug_data.dir == dir_chosen;

h=figure();

ax1 = subplot(2,1,1);

boxchart(ax1, double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title('AUG --> direction = 1 8:00-8:30 segment id vs speed boxplot')

xlabel('segment id kızılay -> polatlı')

ylabel('speed km/h')

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

dir_chosen = 2;

index = aug_data.dir == dir_chosen;

ax1 = subplot(2,1,2);

boxchart(ax1,double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

min_x = min(double(aug_data.segment_id(index)));

max_x = max(double(aug_data.segment_id(index)));

ax1.XLim = [ min_x - 1 , max_x + 1];

ax1.XTick = int64(linspace(min_x,max_x,10));

ax1.XDir = "reverse";

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

xlabel('segment id kızılay -> polatlı')

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title(' AUG <-- direction = 2 8:00-8:30 segment id vs speed boxplot')

time_1 = datetime({'2019-11-18 08:00:00'});

time_2 = datetime({'2019-11-18 10:00:00'});

aux_time = time_1:minutes(30):time_2;

dir_chosen = 1;

h=figure();

for i=1:1:length(aux_time)-1

ax1 = subplot(length(aux_time)-1,1,i);

index = aug_data.dir == dir_chosen & aug_data.time >= aux_time(i) & aug_data.time < aux_time(i+1);

boxchart(ax1,double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

title(ax1,['AUG ', datestr(aux_time(i)), '---', datestr(aux_time(i+1)),' ','direction = ', num2str(dir_chosen), 'segment id vs speed boxplot ndata = ',num2str(sum(index))])

ylim(ax1,[0 200])

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

end

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), -0.1*repmat(yLimits(1),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

set(h,'Units','normalized','Position',[0 0 1 1]);

time_1 = datetime({'2019-11-18 08:00:00'});

time_2 = datetime({'2019-11-18 10:00:00'});

aux_time = time_1:minutes(30):time_2;

dir_chosen = 2;

h=figure();

for i=1:1:length(aux_time)-1

ax1 = subplot(length(aux_time)-1,1,i);

index = aug_data.dir == dir_chosen & aug_data.time >= aux_time(i) & aug_data.time < aux_time(i+1);

boxchart(ax1, double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

ax1.XLim = [ min_x - 1 , max_x + 1];

ax1.XTick = int64(linspace(min_x,max_x,10));

ax1.XDir = "reverse";

title('AUG ', [datestr(aux_time(i)), '---', datestr(aux_time(i+1)),' ','direction = ', num2str(dir_chosen), 'segment id vs speed boxplot ndata = ',num2str(sum(index))])

ylim([0 200])

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

end

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), -0.1*repmat(yLimits(1),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

set(h,'Units','normalized','Position',[0 0 1 1]);

xlabel(' segment id kızılay -> polatlı')

h = figure();

scatter(aug_data.space_diff,aug_data.speed,10,'filled')

xlabel('space difference b/w points (meter)')

ylabel('corresponding speed km/hr')

title(' AUG delta-x vs speed for both dir-s')

h = figure();

scatter(hours(aug_data.time_diff) * 60 * 60,aug_data.speed,10,'filled')

xlabel('time difference b/w points - seconds-')

ylabel('corresponding speed km/hr')

title(' AUG delta-x vs speed for both dir-s')

h = figure();

scatter(aug_data(aug_data.filter,:).space_diff,aug_data(aug_data.filter,:).speed,10,'filled')

xlabel('space difference b/w points (meter)')

ylabel('corresponding speed km/hr')

title({'AUG* delta-x vs speed for both dir-s', ...

' not: uzun yolda segment lengthler uniform gidiyor sanırım', ...

'segment length kısa -> şehiriçi -> hız düşük'})

h = figure();

scatter(hours(aug_data(aug_data.filter,:).time_diff) * 60 * 60,aug_data(aug_data.filter,:).speed,10,'filled')

xlabel(' time difference b/w points - seconds-')

ylabel('corresponding speed km/hr')

title({'AUG* delta-x vs speed for both dir-s', ...

'not: 1/n fonksiyonu v = dx/dt yine de biraz mantıksız? dx sabit mi'})

h = figure();

scatter(hours(aug_data(aug_data.filter,:).time_diff) * 60 * 60,aug_data(aug_data.filter,:).space_diff,10,'filled')

xlabel(' time difference b/w points - seconds-')

ylabel('delta-x')

title({'AUG* delta-x vs delta-t'})

h=figure();

ax1 = subplot(2,1,1);

dir_chosen = 1;

index = (aug_data.dir == dir_chosen) & aug_data.filter;

boxchart(ax1, double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

title(ax1,'AUG* --> direction = 1 8:00-10:00 segment id vs speed boxplot after outlier')

xlabel(ax1,'segment id')

ylabel(ax1,'speed km/h')

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

dir_chosen = 2;

index = (aug_data.dir == dir_chosen) & aug_data.filter;

ax1 = subplot(2,1,2);

boxchart(ax1, double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

min_x = min(double(aug_data.segment_id(index)));

max_x = max(double(aug_data.segment_id(index)));

ax1.XLim = [ min_x - 1 , max_x + 1];

ax1.XTick = int64(linspace(min_x,max_x,10));

ax1.XDir = "reverse";

set(h,'Units','normalized','Position',[0 0 1 0.5]);

title(ax1,'AUG* <-- direction = 2 8:00-10:00 segment id vs speed boxplot after outlier');

xlabel(ax1,'segment id kızılay ->polatlı');

ylabel(ax1,'speed km/h');

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), 0.7*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

ax1.XLim = [min(double(aug_data.segment_id(index))) , max(double(aug_data.segment_id(index)))];

time_1 = datetime({'2019-11-18 08:00:00'});

time_2 = datetime({'2019-11-18 10:00:00'});

aux_time = time_1:minutes(30):time_2;

dir_chosen = 1;

h=figure();

for i=1:1:length(aux_time)-1

ax1 = subplot(length(aux_time)-1,1,i);

index = aug_data.dir == dir_chosen & aug_data.time >= aux_time(i) & aug_data.time < aux_time(i+1) & aug_data.filter;

boxchart(ax1, double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

min_x = min(double(aug_data.segment_id(index)));

max_x = max(double(aug_data.segment_id(index)));

ax1.XLim = [ min_x - 1 , max_x + 1];

% ax1.XTick = int64(linspace(min_x,max_x,10));

% ax1.XDir = "reverse";

title(['AUG* ',datestr(aux_time(i)), '---', datestr(aux_time(i+1)),' ','direction = ', num2str(dir_chosen), 'segment id vs speed boxplot ndata = ',num2str(sum(index))]);

ylim([0 200]);

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

end

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), -0.6*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

set(h,'Units','normalized','Position',[0 0 1 1]);

xlabel('segment id kızılay ->polatlı')

time_1 = datetime({'2019-11-18 08:00:00'});

time_2 = datetime({'2019-11-18 10:00:00'});

aux_time = time_1:minutes(30):time_2;

dir_chosen = 2;

h=figure();

for i=1:1:length(aux_time)-1

ax1 = subplot(length(aux_time)-1,1,i);

index = aug_data.dir == dir_chosen & aug_data.time >= aux_time(i) & aug_data.time < aux_time(i+1) & aug_data.filter;

boxchart(ax1, double(aug_data.segment_id(index)) ,aug_data.speed(index),"MarkerStyle",'.','MarkerColor','red')

min_x = min(double(aug_data.segment_id(index)));

max_x = max(double(aug_data.segment_id(index)));

ax1.XLim = [ min_x - 1 , max_x + 1];

ax1.XTick = int64(linspace(min_x,max_x,10));

ax1.XDir = "reverse";

title(['AUG* ',datestr(aux_time(i)), '---', datestr(aux_time(i+1)),' ','direction = ', num2str(dir_chosen), 'segment id vs speed boxplot ndata = ',num2str(sum(index))]);

ylim([0 200]);

yLimits = get(ax1,'YLim');

line(ax1,[markers.segment_id(markers.dir == dir_chosen),markers.segment_id(markers.dir == dir_chosen)] , yLimits,'LineStyle','--','LineWidth',0.1);

end

text(ax1, double(markers.segment_id(markers.dir == dir_chosen)), -0.6*repmat(yLimits(2),sum((markers.dir == dir_chosen)),1),(markers.names(markers.dir == dir_chosen)),'Rotation',90);

set(h,'Units','normalized','Position',[0 0 1 1]);

xlabel('segment id kızılay ->polatlı')

h = figure();

ax1 = subplot(1,2,1);

dir_chosen = 1;

index = aug_data.filter & aug_data.dir == dir_chosen & ~isnan(aug_data.speed);

color_1 = zeros(sum(index),1);

a = findgroups(aug_data(index,{'arac_id','route_id'}));

% split apply?

for i=1:1:max(a)

aux_index = a == i;

[mat,color_1(aux_index),limits] = get_speed_color(aug_data(index,:).speed(aux_index));

end

scatter(ax1,aug_data.time(index), aug_data.distance_from_start(index),10,mat(color_1,:),'filled');

xlabel(ax1,'time')

ylabel(ax1,'distance from start point POLATLI (m)');

colormap(mat);

caxis(ax1,[0,limits(1)+15]);

c = colorbar;

c.Ticks = linspace(0,limits(1)+15,size(limits,2)+2);

title(ax1,'AUG* dir = 1')

xLimits = get(ax1,'XLim');

line(ax1,[xLimits(1),xLimits(2)],[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] ,'LineStyle','--','LineWidth',0.1,'Color','magenta');

text(ax1, repmat(xLimits(1)- minutes(7),sum((markers.dir == dir_chosen)),1), double(markers.distance_from_start(markers.dir == dir_chosen)),(markers.names(markers.dir == dir_chosen)),'HorizontalAlignment', 'right');

ax1 = subplot(1,2,2);

index = osrm_data_sp.filter & osrm_data_sp.dir == dir_chosen & ~isnan(osrm_data_sp.speed);

color_1 = zeros(sum(index),1);

a = findgroups(osrm_data_sp(index,{'arac_id','route_id'}));

% split apply?

for i=1:1:max(a)

aux_index = a == i;

[mat,color_1(aux_index),limits] = get_speed_color(osrm_data_sp(index,:).speed(aux_index));

end

scatter(ax1,osrm_data_sp.time(index), osrm_data_sp.distance_from_start(index),10,mat(color_1,:),'filled');

xlabel(ax1,'time')

ylabel(ax1,'distance from start point POLATLI (m)');

colormap(mat);

caxis(ax1,[0,limits(1)+15]);

c = colorbar;

c.Ticks = linspace(0,limits(1)+15,size(limits,2)+2);

title(ax1,'MAT* dir = 1')

xLimits = get(ax1,'XLim');

line(ax1,[xLimits(1),xLimits(2)],[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] ,'LineStyle','--','LineWidth',0.1,'Color','magenta');

text(ax1, repmat(xLimits(1)- minutes(7),sum((markers.dir == dir_chosen)),1), double(markers.distance_from_start(markers.dir == dir_chosen)),(markers.names(markers.dir == dir_chosen)),'HorizontalAlignment', 'right');

set(h,'Units','normalized','Position',[0 0 1 0.5]);

h = figure();

ax1 = subplot(1,2,1);

dir_chosen = 2;

index = aug_data.filter & aug_data.dir == dir_chosen & ~isnan(aug_data.speed);

color_1 = zeros(sum(index),1);

a = findgroups(aug_data(index,{'arac_id','route_id'}));

% split apply?

for i=1:1:max(a)

aux_index = a == i;

[mat,color_1(aux_index),limits] = get_speed_color(aug_data(index,:).speed(aux_index));

end

scatter(ax1,aug_data.time(index), aug_data.distance_from_start(index),10,mat(color_1,:),'filled');

xlabel(ax1,'time')

ylabel(ax1,'distance from start point POLATLI (m)');

colormap(mat);

caxis(ax1,[0,limits(1)+15]);

c = colorbar;

c.Ticks = linspace(0,limits(1)+15,size(limits,2)+2);

title(ax1,'AUG* dir = 2')

xLimits = get(ax1,'XLim');

line(ax1,[xLimits(1),xLimits(2)],[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] ,'LineStyle','--','LineWidth',0.1,'Color','magenta');

text(ax1, repmat(xLimits(1)- minutes(7),sum((markers.dir == dir_chosen)),1), double(markers.distance_from_start(markers.dir == dir_chosen)),(markers.names(markers.dir == dir_chosen)),'HorizontalAlignment', 'right');

ax1 = subplot(1,2,2);

index = osrm_data_sp.filter & osrm_data_sp.dir == dir_chosen & ~isnan(osrm_data_sp.speed);

color_1 = zeros(sum(index),1);

a = findgroups(osrm_data_sp(index,{'arac_id','route_id'}));

% split apply?

for i=1:1:max(a)

aux_index = a == i;

[mat,color_1(aux_index),limits] = get_speed_color(osrm_data_sp(index,:).speed(aux_index));

end

scatter(ax1,osrm_data_sp.time(index), osrm_data_sp.distance_from_start(index),10,mat(color_1,:),'filled');

xlabel(ax1,'time')

ylabel(ax1,'distance from start point POLATLI (m)');

colormap(mat);

caxis(ax1,[0,limits(1)+15]);

c = colorbar;

c.Ticks = linspace(0,limits(1)+15,size(limits,2)+2);

title(ax1,'MAT* dir = 2')

xLimits = get(ax1,'XLim');

line(ax1,[xLimits(1),xLimits(2)],[markers.distance_from_start(markers.dir == dir_chosen),markers.distance_from_start(markers.dir == dir_chosen)] ,'LineStyle','--','LineWidth',0.1,'Color','magenta');

text(ax1, repmat(xLimits(1)- minutes(7),sum((markers.dir == dir_chosen)),1), double(markers.distance_from_start(markers.dir == dir_chosen)),(markers.names(markers.dir == dir_chosen)),'HorizontalAlignment', 'right');

set(h,'Units','normalized','Position',[0 0 1 0.5]);

time_1 = datetime({'2019-11-18 08:00:00'});

time_2 = datetime({'2019-11-18 10:00:00'});

aux_time = time_1:minutes(15):time_2;

selected_segment_id = 43;

h=figure();

aux_data = aug_data(aug_data.filter & aug_data.segment_id == selected_segment_id,:);

aux_data_2 = osrm_data_sp(osrm_data_sp.filter & osrm_data_sp.segment_id == selected_segment_id,:);

min_x = min([aux_data.speed;aux_data_2.speed]);

max_x = max([aux_data.speed;aux_data_2.speed]);

for i=1:1:length(aux_time)-1

ax1 = subplot(2,length(aux_time)-1,i);

index = aux_data.time >= aux_time(i) & aux_data.time < aux_time(i+1);

histogram(ax1, aux_data.speed(index),20)

% xlabel(ax1,'bins')

ylabel(ax1,'speed freq')

xlim(ax1,[min_x - 5 , max_x + 5])

ylim(ax1,[0 28])

title(ax1,{datestr(aux_time(i)), ...

datestr(aux_time(i+1)),'AUG*',[ 'segment - id = ',num2str(selected_segment_id)]});

ax2 = subplot(2,length(aux_time)-1,length(aux_time)-1+i);

index = aux_data_2.time >= aux_time(i) & aux_data_2.time < aux_time(i+1);

histogram(ax2, aux_data_2.speed(index),20)

% xlabel(ax2,'')

ylabel(ax2,'speed freq')

xlim(ax2,[min_x - 5 , max_x + 5])

ylim(ax2,[0 28])

title(ax2,'MAT*')

end

set(h,'Units','normalized','Position',[0 0 1 0.5]);

h=figure();

selected_segment_id = 40;

fig_num = 8;

for i=0:1:(fig_num-1)

ax1 = subplot(1,fig_num,i+1);

histogram(ax1, aug_data(aug_data.filter & aug_data.segment_id == selected_segment_id+i,:).speed,20,'FaceAlpha',0.3,'EdgeAlpha',0.7)

hold(ax1,'on');

histogram(ax1, osrm_data_sp(osrm_data_sp.filter & osrm_data_sp.segment_id == selected_segment_id+i,:).speed,20,'FaceAlpha',0.3,'EdgeAlpha',0.7)

% min_x = min(double(aux_data.segment_id(index)));

% max_x = max(double(aux_data.segment_id(index)));

% ax1.XLim = [ min_x - 1 , max_x + 1];

% ax1.XTick = int64(linspace(min_x,max_x,10));

% ax1.XDir = "reverse";

% title([datestr(aux_time(i)), '---', datestr(aux_time(i+1)),' ','direction = ', num2str(aux_dir), 'segment id vs speed boxplot ndata = ',num2str(sum(index))]);

title(['segment -id= ', num2str(selected_segment_id+i),' for all time'])

legend('AUG*','MAT*')

% xlim([10 110]);

end

set(h,'Units','normalized','Position',[0 0 1 0.5]);